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Abstract. Using numerical simulations we study the slow dynamics of a colloidal 
hard-sphere fluid adsorbed in a matrix of disordered hard-sphere obstacles. We 
calculate separately the contributions to the single-particle dynamic correlation 
functions due to free and trapped particles. The separation is based on a Delaunay 
tessellation to partition the space accessible to the centres of fluid particles into 
percolating and disconnected voids. We find that the trapping of particles into 
disconnected voids of the matrix is responsible for the appearance of a nonzero long- 
time plateau in the single-particle intermediate scattering functions of the full fluid. 
The subdiffusive exponent z, obtained from the logarithmic derivative of the mean- 
squared displacement, is observed to be essentially unaffected by the motion of trapped 
particles: close to the percolation transition, we determined z ~ 0.5 for both the full 
fluid and the particles moving in the percolating void. Notably, the same value of 
z is found in single-file diffusion and is also predicted by mode-coupling theory along 
the diffusion-localisation line. We also reveal subtle effects of dynamic heterogeneity in 
both the free and the trapped component of the fluid particles, and discuss microscopic 
mechanisms that contribute to this phenomenon. 
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1. Introduction 

As colloidal fluid particles move within the external, random field induced by a quenched, 
disordered arrangement of particles ("matrix"), they can be either hindered in their 
motion by so-called cages formed by other fluid (mobile) particles, or trapped by the 
immobile particles of the matrix. In the former case, these cages eventually open due 
to the steady, collective motion of the fluid, while in the latter case whether or not a 
fluid particle can escape a certain region of space depends only on pre-set geometric 
and/or energetic restrictions. Each of these processes alone as well as their interplay are 
responsible for slowing down the dynamics of the fluid, leading thereby to a complex 
dynamic behaviour of the system. 

These observations represent yet another confirmation that the properties of a 
fluid that is exposed to a random porous matrix differ substantially from those of 
a bulk fluid. While these phenomena have been extensively studied in experiment 
for a wide variety of systems (e.g., [U El El EH El [6] ), related activities in theory and 
computer simulations considerably lag behind. This discrepancy is mainly due to two 
problems: (i) a realistic parametrisation of typical matrices studied in experiment is 
very difficult, and (ii) the evaluation of the physical properties of the fluid particles is 
conceptually difficult and numerically expensive. As a consequence, so far most, of the 
theoretical and simulation-based studies have focused on static structural [3 El El HQ] , 
and thermodynamic properties [HI [121 H31 HU EEl EES]. Only during the past few 
years, investigations on the dynamic properties of such systems have come within reach 

na nn nn in]. 

Recently, a major breakthrough has been achieved in the realm of theory: 
using concepts from the statistical description of so-called "quenched-annealed" (QA) 
mixtures (see below), a framework has been put forward that allows to determine 
dynamic correlation functions of fluids confined in a porous matrix [201 EH E2]- The 
theory represents an extension to mode-coupling theory (MCT) [23], a successful 
framework that can be used to obtain the dynamic correlation functions of a system 
solely from knowledge about its static structure functions. We shall henceforth refer to 
the theory as "ROZ-MCT", owing to the fact that it is based on mathematical methods 
used in deriving the replica Ornstein-Zernike (ROZ) equations [21], and that it can use 
structural information provided by the ROZ equations as a source of input. 

Meanwhile, thanks to the considerably increased computational power of present- 
day computers, it is possible to systematically investigate the slow dynamic properties 
of fluids confined in a disordered environment in computer models. Simulating such 
systems is costly for multiple reasons: (i) long simulations are necessary to reveal 
the dynamic slowing-down of dense systems (ii) a large body of system parameter 
combinations needs to be considered (iii) many independent realisations of the same 
system have to be generated. The latter requirement arises from the prescription to 
evaluate the physical properties of QA systems, which involves an explicit double average 
procedure: the first average is the usual thermodynamic one to be taken over all possible 
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configurations of fluid particles given a fixed matrix configuration, the second average 
is a disorder average to be performed over all possible matrix representations fulfilling 
pre-imposed criteria. 

To date, the ROZ-MCT framework has been applied only to QA models in which 
hard-sphere particles are confined in a hard-sphere matrix — either with fluid and matrix 
being of equal size [201 [2H [22] or different sizes (25] . Within the QA model, a matrix 
configuration is obtained by taking a snapshot of an equilibrated fluid (quenched 
component); subsequently the particles of the fluid (annealed component) are immersed 
in this configuration. Shortly after publication of the theoretical predictions, extensive 
and systematic molecular dynamics (MD) simulations have been performed for precisely 
this system [25 1 |26 | [27]. While simulations confirmed most of the intriguingly complex 
features that the theory predicted for the kinetic diagram (such as a discontinuous and 
a continuous glass transition as well as a continuous diffusion- localisation transition), 
others are yet to be found (in particular the re-entrant glass transition scenario). In this 
context, it is worth mentioning other recent works which investigate the dynamics of a 
fluid [HI [191 [28] or of tracer particles [291 EH Ell E2] moving in a matrix of quenched 
particles. 

Our previous work [251 I2S] was dedicated to an "overall" investigation of the 
dynamic properties (in terms of both single-particle and collective dynamic properties) 
of hard spheres confined in a disordered matrix of hard-sphere particles, and to a 
comparison of these results to the theoretical predictions [201 [2T| [22]. Computer 
simulations offer the distinct possibility to proceed one step further: in contrast to 
the theoretical framework, in simulations it is possible to distinguish between caging 
and trapping by geometrically analysing the voids formed by the matrix particles. 

The pertinent procedure reads as follows. Given a particular realisation of a 
state point, the first step is to classify the voids formed by the matrix particles while 
disregarding the fluid particles. The algorithm used is detailed in Ref. [33]; in short, we 
distinguish between disconnected voids of finite volume from which — due to geometric 
restrictions — particles cannot escape, and the percolating void that — taking into account 
the periodic boundary conditions — extends over the entire space and in which particles 
can propagate infinitely far. (Note that multiple percolating voids can only exist in 
systems of finite size.) In the second step, for each fluid particle in the system realisation 
we use this void analysis to determine in which class of void the particle is located. The 
fluid particles thus fall in two classes: "trapped" particles (superscript "trap") that 
populate the disconnected voids, and "free" particles (superscript "free") that move 
in the percolating void. We then compute a set of single-particle dynamic correlation 
functions [namely the self-intermediate scattering function F s (k,t), the mean squared 
displacement 6r 2 (t), and the self-part of the direction- averaged van Hove function 
G s (r,t)} separately for the trapped and the free particles [i.e., F s tmp (/c,t) and F s free (/c, t), 
etc.]. 

In an effort to specify and to quantify the respective contributions of trapped and 
free particles to the dynamic structure, F s (k, t) and Sr 2 (t) are discussed along two specific 
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pathways in the system's parameter space. The latter is spanned by the packing fraction 
of the fluid, 0f, and that of the matrix, </> m . The properties of G s (r,t) for the trapped 
and the free particles are discussed for selected state points in the (0 m , <f>t) plane. 

Notwithstanding the presentation of quantitative results in Sec. [31 it is worth 
clarifying what qualitative behaviour should be expected for the dynamic correlation 
functions. In the long-time limit, the contributions of the trapped and the free particles 
differ markedly, which bears physical significance. While the single-particle intermediate 
scattering function of the free particles will always relax completely for t — > oo, this is 
never the case for F s trap (fc, t). On the other hand, the mean-squared displacement of the 
free particles will always recover diffusive behaviour while Srf (t) will always saturate 
as t — >• oo, with the limiting value reflecting the average pore size. Finally, for the 
free particles the self-part of the van Hove function will rapidly relax towards zero for 
distances in the vicinity of zero. Conversely, Gg rap (r~0,t) will not decay below a value 
significantly above zero. Instead, it will approach a steady-state distribution that reflects 
the distribution in trap sizes and shapes. 

The paper is organised as follows: in the subsequent section (Sec. [2]) we briefly 
summarise the model and our simulation technique (referencing previous publications 
for more details). Section [3] is dedicated to the results: after a brief summary of our void 
identification algorithm as well as a short discussion of the void analysis of the state 
points treated in this contribution, we present in detail our results for the single-particle 
correlation functions. In Sec. 0] we discuss implications of those results, and present 
explanations for the observations. The paper is closed with concluding remarks. 

2. Model and theoretical approaches 

Following our previous work, we have used the QA protocol to describe the properties 
of fluid particles that are exposed to an external field of a random porous matrix. In 
this picture, the fluid particles move within the space left free by the matrix: both types 
of particles interact with each other, but the mobile fluid particles cannot displace the 
matrix particles. Choosing all particles to interact as hard spheres, any state point is 
fully characterised by a combination of <ft m and 0f. Different, but equivalent matrix 
configurations of a given state point are obtained from an equilibrated hard-sphere fluid 
at a prescribed m . For this, a pertinent simulation is halted at different times and the 
particle positions are subsequently used as positions of the matrix particles in a QA 
system. In all MD simulations we have used an event-driven MD algorithm adapted to 
keeping the matrix particles at fixed positions. Periodic boundary conditions and the 
minimum image convention were employed. 

As prescribed by thermodynamics [3U ESI El], observables are obtained from a 
double- averaging procedure: the first trace, taken over the degrees of freedom of the 
fluid particles for a particular matrix configuration, is realised using a time average 
along the simulation run; the second one, taken over different yet equivalent matrix 
configurations, is realised as an ensemble average over independent simulation runs. 
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All state points were studied via ensembles in which the number of fluid particles (TVf) 
and the number of matrix particles (N m ) added up to at least 1000. For systems with 
elevated matrix packing fractions this number was considerably enhanced in an effort 
to guarantee that for all state points considered the ensemble contain at least 50 fluid 
particles. A method for finding an initial configuration of QA systems at large values of 
the total packing fraction, </> m + (f> { , has been presented and discussed in the Appendix 
of Ref. [26]. 

In the first phase of the simulation the system is equilibrated; for clarifying 
comments on this delicate issue we refer the reader to Section II of Ref. [26]. In 
case the system was successfully brought to equilibrium, we subsequently recorded 
the particle positions in a production run extending over the same simulation time 
as the equilibration run. We typically allowed for a maximum of 30 000 r for the 
equilibration run; for selected systems this figure was increased to the tenfold value. Here 
t = 'ma 2 1 'k^T is the unit of time [36], where the particle mass m, the temperature T, 
and the particle diameter a were chosen to be unity. 

The definitions of all observed quantities presented in this work [F s (k,t), Sr 2 (t), 
z(t), and G s (r,t)] have been compiled in Ref. [26], to which we refer the reader. 

3. Results 

3.1. Void analysis 

Due to the statistic nature of our matrix model, at any value of the matrix packing 
fraction <fr m there is a nonvanishing probability for the immobile matrix particles to form 
traps - void spaces from which fluid particles enclosed therein are not able to escape due 
to geometric restrictions. At low values of m the number of traps is small, and most 
fluid particles can move infinitely far away from their initial location since the void in 
which they are located extends (taking into account the periodic boundary conditions) 
over the entire space. The latter volume, if present, is termed the "percolating" void. At 
sufficiently high m , on the other hand, all voids will be of finite size and represent traps 
from which the fluid particles cannot escape. The transition between the two scenarios 
occurs at the so-called percolation transition. (For an overview over the large field 
of percolation transitions see, e.g., [37]). For an infinitely-large system this transition 
is sharp, i.e., it occurs at some well-defined value 0^ of the matrix packing fraction. 
However, if the system is represented by a finite number of particles then this transition 
is smeared out; for this reason in computer simulations matrix realisations featuring a 
percolating void are encountered also beyond 0^. 

We employed an algorithm based on a Delaunay tessellation — details are presented 
in Ref. [33] — in order to identify for a given matrix configuration whether a void 
constitutes a trap or the percolating void. Using the position of a fluid particle at 
an arbitrary instance of time, it is thus possible to determine whether the particle in 
question is "free" or "trapped" (according to the previously-introduced notion). The 
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Figure 1. Visualisation of the procedure to identify trapped and free fluid particles, 
(a) Dark grey spheres: matrix particles. Volumes in lighter (red) shade: voids (volumes 
accessible to the fluid particles), (b) Same as previous, with the matrix particles 
removed, (c) Same as previous, with the percolating void highlighted in a light (green) 
shade, (d) Same as previous, with a sample set of fluid particles inserted into the voids. 
Light (grey) spheres: free particles. Dark (blue) spheres: trapped particles. 
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Figure 2. Parameter space for QA systems of hard spheres, spanned by the packing 
fraction of the matrix particles, </> m , and that of the fluid particles, <f>{. The paths 
and points relevant to this contribution are indicated as large arrows and crosses. 
For reference, they are superimposed to data taken from Fig. 4 in Ref. [55], in which 
simulations were performed at state points indicated by the solid circles, and systems 
right of the dashed line were defined to be "arrested" based on F s (k=7, t). 



method, which has recently received considerable attention in contexts similar to the 
one of this work j32j [381 EH] 5 is sketched in Fig. CD 

3.2. Dynamic correlation functions 

Making use of the distinction in trapped and free particles, the single-particle correlation 
functions F s (k,t), 5r 2 (t), and G s (r,t) can be split up in an additive way into two 
component-weighted contributions: one stemming from the iV f trap trapped particles 
(index "trap"), the other one from the iV f free free particles (index "free"). Since the 
fraction iV f trap /(iV f free + iV f trap ) increases with m , the correlation functions of the trapped 
particles will suffer from poor statistics low m , whereas for <p m > (p* m the results for the 
correlation functions of the free particles will be affected by a relatively large numerical 
error. 

To facilitate comparison with previous data j25j [26] , in this contribution we discuss 
the results in a similar fashion. This includes the wave vector at which we examine 
F s (k,t), namely k = 7 [this value is close to the first peak in the structure factor S(k)], 
as well as the choice of specific paths in the parameter space — depicted in Fig. [2] — at 
which we evaluated F s (k,t) and Sr 2 (t). Path I in these works (extending at fixed, low 
m while increasing 0f) will not be considered here since trapping — the key property 
studied in this contribution — plays only a minor role in this region of the kinetic diagram. 




Figure 3. Single-particle intermediate scattering function as a function of time t at 
fixed k = 7 and <f> m = 0.20 for a series of </>f. (a) for the free particles, F* lee (k,t), (b) 
for the trapped particles, F^ ra,p (k,t). Error bars represent one standard deviation of 
the mean for different system realisations. 



Path II denotes systems at the intermediate fixed value m = 0.2 and increasing fluid 
packing fractions; note that the matrix packing fraction along path II is lower than </>^ 
[251 1261 [33] . Finally, along path III <p{ is kept fixed at 0.1 while increasing 0f over an 
interval that includes the percolation transition of the voids at 0^. 

We start with the discussion of F s (k,t) along path II; the respective contributions 
of the free and the trapped particles are displayed in Fig. El Essentially, F s frcc (fc=7, t) 
exhibits the same behaviour as the total correlator, F s (k=7,t) (shown in Fig. 8(b) 
of [26]), except for the fact that the long-time tail now tends towards zero for all 
equilibrated samples. In contrast, F s trap (/c=7, t) shows a long-time plateau that increases 
with 0f, which can be explained as follows: since we keep the value of m fixed, the void 
size distribution is the same for all systems along path II. The increase in the plateau 
is therefore due to the increasing number of trapped fluid particles, which cause mutual 
jamming effects inside the trapping voids that in turn prevent the correlator from fully 
relaxing as t tends to infinity. 

The consequences of an increasing fluid packing fraction are even more obvious 
in the respective contributions to the mean-squared displacement, 5r 2 (t), displayed 
in Fig. HI In the long-time limit, Srf rec (t) shows normal diffusive behaviour for all 
equilibrated systems. In contrast, in the long-time limit of 5rf rap (t) a saturation is 
observed; the plateau values decrease as 0f is increased, which is due to the crowding 
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Figure 4. Mean-squared displacement as a function of time t at fixed fa n = 0.20 for a 
series of fa. (a) for the free particles, <frf ree (i), (b) for the trapped particles, <5rf rap (£). 
Error bars: see Fig. |3] 




Figure 5. Logarithmic derivative of the mean-squared displacement as a function of 
time t at fixed fan = 0.20 for a series of fa. (a) for the free particles, Zf iee (t), (b) for 
the trapped particles, z trap (£). 




Figure 6. Single-particle intermediate scattering function as a function of time t at 
fixed k = 7 and fa = 0.10 for a series of </> m . (a) for the free particles, F s free (fc, i), (b) 
for the trapped particles, F s trap (fc, £). Error bars: see Fig. [3] 



to which the fluid particles are subjected for larger fluid packing fractions. Assuming 
an instantaneous power-law behaviour 5r 2 (t) oc t z , the respective long-time behaviour 
of both contributions to the mean-squared displacement is nicely reflected in the 
subdiffusion exponent z(t), as displayed in Fig. |5j it quickly drops to zero for the 
trapped particles, whereas for the free particles Zf Tee (t) tends to unity as t — > oo, which 
corresponds to diffusive behaviour (for the three lowermost curves in Fig. 0(a) the limit is 
not approached yet). Notably, for several decades in time an intermediate-time plateau 
can be observed in Zf ree (t); the value of this plateau is virtually the same as that found 
for the exponent pertaining to 5r 2 (t) for the full fluid (see Fig. 9(b) in [26]). 

The contributions to F s (k,t) along path III are displayed in Fig. El Disregarding in 
panel (a) the state point above the percolation threshold (which is poorly equilibrated 
with respect to the free particles), we find that F^ ree (k=7,t) relaxes all the way to zero 
for t — > oo. This contrasts with the single-particle intermediate scattering function of 
the full fluid — cf. Fig. 10(b) in Rcf. [26J — which attains a nonzero long-time plateau 
that increases with the matrix packing fraction. Hence, the nondecaying part of F s (k, t) 
must be entirely due to the trapped particles; this fact is reflected in the significant 
nonrelaxing part of F s trap (/c, t). We observe that _F s trap (/c, t— >oo) exhibits a nonmonotonic 
behaviour, attaining a minimum value at <p m ~ 0.25, which is close to the percolation 
transition in the void space [33]. Sr 2 Tee (t) and Zf Tee (t) — the latter is presented in Fig. [El- 
show features similar to those of their equivalents for the full fluid (cf. Fig. 11 in |26j). 
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Figure 7. Mean-squared displacement as a function of time t at fixed </>f = 0.10 for a 
series of <p m . (a) for the free particles, Srf Tee (t), (b) for the trapped particles, 6rf (i). 
Error bars: see Fig. |3] 
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Figure 8. Logarithmic derivative of the mean-squared displacement as a function of 
time t at fixed <j>f = 0.10 for a series of <f> m . (a) for the free particles, Zf iee (t), (b) for 
the trapped particles, z tI&p (t). 




Figure 9. Self-part of the direction-averaged van Hove function as a function of r at 
point "B", denned via 4> m = 0.20 and fa = 0.22. (a) for the free particles, G l s loe (r, t), 
(b) for the trapped particles, G' rap (r, t). Curves are approximately logarithmically 
spaced in time; some are highlighted for visual guidance. Error bars: see Fig. [3] 



apart from two differences: firstly, Z{ Tee (t) shows an onset to recover diffusive behaviour 
for m = 0.2625, as opposed to z(t) which shows a trend to saturation. Therefore, in 
the time range in question (t ~ 10 5 ) the mean-squared displacement of the full fluid has 
to be dominated by the contribution of the trapped particles. The second interesting 
feature is the much clearer subdiffusive regime for </> m = 0.2625, which now remains at 
an essentially invariable Z{ ree ~ 0.5 for almost three decades in time. The mean-squared 
displacement of the trapped particles reveals a similarly interesting behaviour: while 
<5rf rap (t), shown in Fig. [71 is subdiffusive with z trap (t) < 0.5 at all times, for matrix 
packing fractions close to the percolation threshold Zt ra , p (t) attains much larger values 
of almost 0.5 at intermediate times, and 5rf (t) saturates at much larger values in this 
m regime. The latter indicates that the average trap size attains a maximum at </>^. 
This behaviour is well known from percolation theory, which predicts the average trap 
size to diverge; onsets of this phenomenon were also shown to exist in Refs. [331 132] . 

The self-part of the direction- averaged van Hove function, G s (r, t), shall be discussed 
for two selected state points, labelled point "B" and point "C" , which are indicated in 
Fig. Point B is located on path II (fixed m = 0.20) at a fluid packing fraction 
0f = 0.22, where the latter value has been shown sufficient to cause a relatively slow 
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Figure 10. Two-dimensional schematic of a long dead-end channel. Black solid circles: 
matrix particles that form the dead-end channel. Grey (coloured) solid circles: fluid 
particles in the channel. Light grey hatched circles: other matrix particles. Light 
grey open circles: other fluid particles. From light green to deep blue (colour online) 
there is a decreasing probability for the particle to leave its position via collective 
rearrangements. 



relaxation of the fluid's dynamic features [251 [26] . The respective contributions to the 
van Hove functions are depicted in Fig. [9j The main peak of G^ ee (r, t) at small distances 
drops rapidly at small times. However, for large times it is still present; this is probably 
due to fluid particles that are located at the end of long, narrow dead-end channels like 
the one schematically drawn in Fig. [TDJ Those channels are connected to the percolating 
void, but very improbable collective rearrangements would have to occur for the end- 
particle to leave its position. At intermediate times and r ~ 1.1 we see a pronounced 
peak in G^ cc (r, t), which indicates hopping processes of the free particles. Gg rap (r, t), on 
the other hand, quickly reaches a stationary distribution; the latter is relatively narrow 
which indicates that for m = 0.20 only a small fraction of pores is large enough to 
accommodate more than one particle. For large times the van Hove function displays a 
small shoulder at r ~ 0.5, reflecting the pore size distribution; the small peak at r ~ 1 
emerging for large times evidences that despite the small average trap size there exist 
traps large enough for particles to collectively rearrange. 

Point C, the other point at which we examined G s (r,t), is located on path III 
and is characterised by 0f = 0.1 and <p m = 0.25; notably, the latter figure is close to 
the percolation threshold. The respective contributions to the van Hove functions are 
depicted in Fig. [TTJ The main peak of G^ ee (r,t) at r ~ drops to zero as t tends to 
infinity. For r = a minute peak seems to persist; it might be associated with the 
trapping phenomena discussed in Sec. HI Other than a hopping peak at r ~ 1 (less 
pronounced than that in Fig. [9]) there are no unusual features in Gj; ree (r, t). In contrast, 
for this state point G* rap (r, t) shows a behaviour which is distinctively different from that 
of the same correlator for the other state point: at point C the evolution of the stationary 




Figure 11. Self-part of the direction- averaged van Hove function as a function of r at 
point "C", denned via </> m = 0.25 and 4>f = 0.10. (a) for the free particles, G l s loe (r, t), 
(b) for the trapped particles, G' rap (r, t). Curves are approximately logarithmically 
spaced in time; some are highlighted for visual guidance. Error bars: see Fig. [3] 



state at t — > oo takes considerably more time than at point B (cf. Fig. |9]). Also, the 
pronounced peak at small distances in the limiting function evidences a considerable 
amount of trapping; the broad side peak at r ~ 1.1 reflects hopping processes inside the 
traps. Finally, it is worth noting that G7* rap (r, t— >oo) exhibits a very long tail in r; this 
is another indication of the fact that close to the percolation transition there exist very 
large traps. 

4. Discussion 

Our analysis of the dynamic correlation functions for the free and the trapped particles 
provides some deeper insight into the complex relaxation of a colloidal fluid adsorbed in a 
dense immobile matrix. Comparison with results for the full fluid [251 [26] strengthens our 
previous interpretation of the dynamic data and allows to assess the various predictions 
of MCT [201 El [22] from a broader perspective: 

(i) The complex relaxation of the single-particle correlators along a pathway 
representative of the behaviour at intermediate matrix and large fluid densities (path II) 
arises from a genuine superposition of trapping effects and a collective caging mechanism: 
the former leads to a finite long time plateau in F s (k,t), while the latter is responsible 
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for the inflection in the intermediate time regime of this correlation function. This 
observation strongly supports our previous interpretation for the data of the the full 
fluid along this path [251 EE] and indicates that this behaviour may correspond, at least 
on a qualitative level, to the crossing of the two arrest lines (diffusion-localisation and 
glass-transition) as predicted by MCT. 

(ii) The values of the subdiffusion exponent z for the full fluid [251 EE] clearly show 
that the free fluid particles undergo a subdiffusive process on their own right. Moreover, 
the values of z for the full fluid are only marginally affected by the contribution of the 
trapped particles: our calculation of Zf rec at a packing fraction close to the estimated 
percolation transition is, within error margin, indistinguishable from the value 0.5 found 
in the full system. As an observational fact, the same value of z is featured in single-file 
diffusion [401 l4~Tj and is also predicted by MCT along the diffusion-localisation transition 
[22] • However, it is well-known that MCT cannot properly describe the diffusion of a 
tracer particle in a disordered frozen matrix such as in the Lorentz gas model [221 H2]- 
In particular, in the Lorentz gas the subdiffusion exponent z = 0.32 is observed at the 
percolation transition [291 130] . which differs from the value z = 0.5 predicted by MCT 
at the diffusion-localisation transition in QA systems. Motivated by this observation 
and by related findings in previous works [23 S3] , we have examined the case of non- 
interacting tracer particles moving in a QA-type matrix at the percolation transition 
[44] . We found that in this case z for the full fluid assumes a value around 0.3, which is 
more similar to the value found in the Lorentz gas than to the value predicted by MCT. 
The apparent disagreement between the case of noninteracting tracers and the case of 
a dense fluid moving in such matrices will be subject to future work [4"4"] . 

(iii) The splitting procedure effectively removes the nonzero long-time limit of 
F s (k,t). This method therefore provides a microscopic basis to the empirical procedure 
described in [26], where an ad- hoc subtraction of the long time plateau was performed. 
The agreement between the results of the two procedures (not discussed here) lends 
support to our discussion of the modified kinetic diagram for single-particle properties 
presented in [26] . 

Even more intriguing is the observation of nontrivial dynamic heterogeneity in both 
the trapped and the free component of the fluid particles. For the free particles, the 
van Hove function displayed in Fig. UTT a) evidences that a non-negligible fraction of the 
particles undergo a very slow relaxation. At least three physical mechanisms explaining 
this finding can be identified, (i) Free particles may be located at the end of long 
dead-end channels as shown in Fig. [10J the longer such a channel is, the less likely 
are the collective rearrangements required for a particle to escape from the channel's 
end. (ii) Free particles may mutually jam in a pore as schematised in Fig. PTWa); such 
particles are rigorously confined to that pore since none of them can reach the channel 
that connects the pore to the percolating void, (iii) A trapped particle may act as an 
"effective" matrix particle [see Fig. [T2T b)]: if the matrix restricts its movement to much 
less than its radius, it can block a pathway in the percolating void; this may render the 
void at one side of that pathway a trap. Visual inspection of animated trajectories of 
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Figure 12. Two-dimensional schematics for mechanisms of effective trapping, induced 
by the presence of multiple fluid particles in the system. Black solid circles: relevant 
matrix particles. Dark purple solid circles: trapped particles. Dark blue solid circles: 
effectively-trapped particles. Green open circles: positions a fluid particle can occupy. 
Red open circles: positions a fluid particle can not occupy due to overlaps. Yellow: 
area constituting the effective trap. Light grey hatched circles: other matrix particles. 
Light grey open circles: other fluid particles populating the percolating void, (a) The 
effectively-trapped particles are located in a pore connected to the percolating void; 
however, due to jamming the particles cannot rearrange such that one of them would 
reach the connecting channel, (b) The effectively-trapped particles are located in a 
pore connected to the percolating void; however, another fluid particle (purple) is 
trapped between matrix particles in a way that it blocks the connecting channel for 
the particles in the pore (blue). 



our systems confirmed that all of these processes are present and relevant in the portion 
of the kinetic diagram corresponding to large 0f and m . Note that effects (ii) and (iii) 
arise from the QA-protocol requirement that all available volume (including traps) be 
populated with fluid particles. 

Notably, dynamic heterogeneity is also present among the trapped particles; they 
have been observed to rearrange collectively for instance by exchanging their positions 
in a closed loop. Thus, for the trapped particles an apparent weak relaxation appears 
in the single-particle correlators and a tail emerges at large distances in their van Hove 
function. A more systematic analysis of these striking dynamic features is currently 
underway 

5. Conclusions 

We have presented a computational study concerned with the dynamic properties 
of a model colloidal fluid adsorbed in a disordered, rigid array of obstacles. The 
original aspect of this study consists in the splitting of dynamic correlators into two 
contributions, one originating from the free and the other from the trapped fluid 
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particles. The notion of "free" and "trapped" refers to the structure of the voids 
formed by the matrix: free particles are located in the infinitely-large percolating void 
whereas trapped particles reside in disconnected voids. We compared the relaxation 
patterns of the two components with the corresponding patterns for the full fluid, which 
allows to disentangle the dynamic effects due of crowding and confinement. Further, 
our procedure provides some deeper insight into the relaxation processes at play at the 
microscopic level. Our analysis confirmed the striking superposition of two relaxation 
mechanisms — trapping and caging — at a sufficiently large density of both the fluid 
and the matrix component, and supported the view that MCT correctly captures the 
qualitative features of the dynamics in the corresponding portion of the kinetic diagram. 
Our numerical procedure also unveils that an unexpected dynamic heterogeneity is 
present in the motion of fluid particles within the matrix structure: on the one hand, 
even particles residing in a percolating void may be effectively trapped; on the other 
hand, residual relaxation can take place within disconnected pores. We believe these 
aspects to be relevant also for the assessment of dynamic heterogeneity in more general 
models of fluids adsorbed in porous media. Whether these complex dynamic phenomena 
are present also in systems with different protocols to generate the matrix [271 H5] . soft 
interactions [25], or systems in which a large dynamic asymmetry replaces the quenched 
disorder [361 H7J H2] is a question that remains to be addressed in future studies. 
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